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REMOTE SENSING OF EARTH TERRAIN 


Under the sponsorship of the NASA Contract NAGW-1617, we published 56 refereed jour- 
nal and conference papers for the research on remote sensing of earth terrain. The following is a 
summary of recent research findings. 


Application of Theoretical Models to Active and Passive Remote Sensing of Saline Ice 

The random medium is used to interpret the polarimetric active and passive measurements 
of saline ice. The ice layer is described as a host ice medium embedded with randomly distributed 
inhomogeneities, and the underlying sea water is considered as a homogeneous half-space, the 
scatterers in the ice layer are modeled with an ellipsoidal correlation function. The orientation of 
the scatterers is vertically aligned and azimuthally random. The strong permittivity fluctuation 
theory is employed to calculate the effective permittivity and distorted Born approximation is used 
to obtain the polarimetric scattering coefficients. We also calculate the thermal emissions base on 
the reciprocity and energy conservation principles. The effects of the random roughness at the 
air-ice, and ice- water interfaces are accounted for by adding the surface scattering to the volume 
scattering return incoherently. The above theoretical model, which has been successfully applied 
to analyze the radar backscatter data of the first-year sea ice near Point Barrow, AK, is used to 
interpret the measurements performed in the CRRELEX program. 

Radiative Transfer Theory for Polarimetric Remote Sensing of Pine Forest 

The radiative transfer theory is applied to interpret polarimetric radar backscatter from 
pine forest with clustered vegetation structures. To take into account the clustered structures with 
the radiative transfer theory, the scattering function of each cluster is calculated by incorporat- 
ing the phase interference of scattered fields from each component, subsequently, the resulting 
phase matrix is used in the radiative transfer equations to evaluate the polarimetric backscattering 
coefficients from random medium layers embedded with vegetation clusters. Upon including the 
multi-scale structures, namely, trunks, primary and secondary branches, as well as needles, we 
interpret find simulate the polarimetric radar responses from pine forest for different frequencies 
and looking angles. The preliminary results are shown to be in good agreement with the measured 
backscattering coefficients at the Landes maritime pine forest the MAESTRO-1 experiment. 


- 1 - 




Scattering of Electromagnetic Waves from a Dense Medium Consisting of Correlated Mie Scatter- 
ed with Size Distributions and Applications to Dry Snow 

The scattering of a plane electromagnetic wave obliquely incident on a layer of dense 
medium consisting of dielectric spherical particles of finite sizes and with size distributions is stud- 
ied. The spherical particles are of sizes comparable to wavelength so that Mie scattering is used 
to describe the single particle scattering characteristics. The coherent wave is studied with qua- 
sicrystalline approximation using the cross pair distribution functions of multiple size governed by 
Percus-Yevick approximation. The incoherent scattered wave is calculated with the distorted Born 
approximation with the result expressed in terms of a product of the T-matrices of particles of 
different sizes and permittivities and the Fourier transform of the cross pair distribution functions. 
The coherent wave effective propagation constants, the attenuation rates and the backscattering 
coefficients are illustrated numerically, with examples chosen to illustrate microwave and millime- 
ter wave scattering from snow cover in the frequency range of 5 GHz to 95 GHz, and mean grain 
radius between 0.03 cm to 0.09 cm. Salient features of the numerical results for scattering from 
snow with size distribution are: (1.) Correlated dense medium scattering is less than independent 
scattering at low frequency, a fact that is consistent with controlled laboratory experiment. (2.) 
Scattering from dense medium of a Rayleigh size distribution with an average radius can be much 
larger than the case of monodisperse particles of the sizes identical to that average radius. (3.) The 
scattering attenuation rate increases rapidly with frequency at low frequency regime and begins to 
level off at high frequency regime. (4.) The coherent wave scattering attenuation rate is large in 
snow at frequencies above 15 GHz indicating large optical thickness and the important of multiple 
scattering. Comparisons with extinction measurements of dry snow at 18 GHz, 35 GHz, 60 GHz, 
and 90 GHz are made. 

Variance of Phase Fluctuations of Waves Propagating through a Random Medium 

As an electromagnetic wave propagates through a random scattering medium, such as a 
forest, its energy is attenuated and random phase fluctuations are induced. The magnitude of the 
random phase fluctuations induced is important in estimating how well a Synthetic Aperture Radar 
(SAR) can image objects within the scattering medium. The two-layer random medium model, 
consisting of a scattering layer between free space and ground, is used to calculate the variance 
of the phase fluctuations induced between a transmitter located above the random medium and 
a receiver located below the random medium. The scattering properties of the random medium 
are characterized by a correlation function of the random permittivity fluctuations. The effective 
permittivity of the random medium is first calculated using the strong fluctuation theory, which 
accounts for large permittivity fluctuations of the scatterers. The distorted Bora approximation is 
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used to calculated the first-order scattered field. A perturbation series for the phase of the received 
field in the Rytov approximation is then introduced and the variance of the phase fluctuations 
is solved to first order in the permittivity fluctuations. The variance of the phase fluctuations is 
also calculated assuming that the transmitter and receiver are in the paraxial limit of the random 
medium, which allows an analytic solution to be obtained. Results are compared using the paraxial 
appromination, scalar Green’s function formulation, and dyadic Green’s function formulation. The 
effects studied are the dependence of the variance of the phase fluctuations on receiver location in 
lossy and lossless regions, medium thickness, correlation length and fractional volume of scatterers, 
depolarization of the incident wave, ground layer permittivity, angle of incidence, and polarization. 

Theoretical Modeling for Passive Microwave Remote Sensing of Earth Terrain 

Historically, the volume scattering effects stimulated the development of the continuous 
random medium model and the random discrete scattering model for the description of the media. 
Theoretical treatments were developed along two different paths. Invoking the principle of reci- 
procity, the wave theory based on Maxwell’s equations has been used to calculate the emissivity. 
The other approach was to start with the radiative transfer equations and solved for the brightness 
temperatures directly. Attempts have been made to derive the radiative transfer theory from the 
wave theory. At the same time, both theoretical approaches have been used to calculate the radio- 
metric emissions and to interpret experimentally measured data. The successful interpretation of 
the Cosmos 243 data was perhaps the first most important step towards a serious development of 
the continuous random medium model to account for the volume scattering effects of snow ice fields. 
Subsequent interpretation of measurement results from snow field with both passive radiometers 
and active radar systems established a unique position for its description of earth terrain media. 
Recent efforts in classifying sea ice with correlation function characterization are demonstration of 
acceptance of this model. Future inverse scattering developments will perhaps rely heavily on this 
model. In this paper, we shall illustrate the development of the theoretical models and present data 
matching results with measurements made in snow ice fields and vegetation canopies. The emis- 
sivity calculations for periodical rough surfaces will also be presented and compared to measured 
data. Recent development in polarimetric active remote sensing with synthetic aperture radar has 
created significant theoretical results and practical applications. In passive remote sensing, the 
third and the fourth Stokes parameters for earth remote sensing have not received much attention 
in the past partly due to the expected small values any measurement can yield. We have made ini- 
tial calculations and experimental measurements to show that at least the third Stokes parameter 
can give appreciable number in both theoretical prediction and actual experimentation in a plowed 
field. Such results may have practical implications in measuring wind directions in ocean waves, for 
instance. We believe polarimetric passive remote sensing is a viable field which should be explored 
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in light of its potential applications that may derived from the full Stokes vector instead of its 
first two parameters, i.e.,ters, i.e., the horizontally and vertically polarized brightness temperature 
components. 

Polar imetric Signatures of a Canopy of Dielectric Cylinders Based on First and Second Order 
Vector Radiative Transfer Theory 

Complete polarimetric signatures of a canopy of dielectric cylinders overlying a homoge- 
neous half space are studied with the first and second order solutions of the vector radiative transfer 
theory. The vector radiative transfer equations contain a general nondiagonal extinction matrix 
and a phase matrix. The energy conservation issue is addressed by calculating the elements of the 
extinction matrix and the elements of the phase matrix in a manner that is consistent with energy 
conservation. Two methods are used. In the first method, the surface fields and the internal fields 
of the dielectric cylinder are calculated by using the fields of an infinite cylinder. The phase matrix 
is calculated and the extinction matrix is calculated by summing the absorption and scattering to 
ensure energy conservation. In the second method, the method of moments is used to calculate the 
elements of the extinction and phase matrices. The Mueller matrix based on the first order and 
second order multiple scattering solutions of the vector radiative transfer equation are calculated. 
Results from the two methods are compared. The vector radiative transfer equations, combined 
with the solution based on method of moments, obey both energy conservation and reciprocity. 
The polarimetric signatures, copolarized and depolarized return, degree of polarization, and phase 
differences are studied as a function of the orientation, sizes, and dielectric properties of the cylin- 
ders. It is shown that second order scattering is generally important for vegetation canopy at C 
band and can be important at L band for some cases. 

Branching model for vegetation 

A branching model has been proposed for the remote sensing of vegetation. The frequency 
and angular response of a two-scale cylinder cluster are calculated to demonstrate the significance 
of vegetation architecture. The results indicate that it is necessary for theoretical remote sensing 
models to take into account the architecture of vegetation which plays an important role in deter- 
mining the observed coherent effects. A two- scale branching model is implemented for soybean with 
its interned structure and the resulting clustering effects considered. Furthermore, at the scale of 
soybean fields, the relative locations of soybean plants is described by a pair distribution function. 
The polarimetric backscattering coefficients are obtained in terms of the scattering properties of 
soybean plants and their pair distribution. Theoretical backscattering coefficients are evaluated 
using two pair- distribution functions: independent scatterer and hole correction. Backscattering 
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coefficients calculated by using independent-scatterer pair distribution are in good agreement with 
extensive data collected from soybean fields except for the data near normal incidence. It is found 
that the hole- correction approximation, which prevents two soybean plants from overlaying each 
other, is more realistic. By introducing destructive interference at small angles of incidence, it 
greatly improves the agreement between the model and these three data sets near normal inci- 
dence. Extension to a multi-scale branching model can be achieved by recursion of the two-scale 
modelling approach reported in this research. 

In the past, when radiative transfer theory was applied to the modeling of vegetation, the 
average phase matrix of vegetation layer was approximated by an incoherent sum of the phase 
matrices of individual vegetation elements. However, it can be observed that most of vegetation 
species exhibits a branching structure, whose significance has been demonstrated by the frequency 
and angular response of a two- scale cylinder cluster. This indicates that it is necessary for theoretical 
remote sensing models to account for the architecture of vegetation. 

In this research, radiative transfer theory is applied to vegetation with clustered structures. 
To take into account vegetation structure in the radiative transfer theory, the phase matrix of a 
vegetation cluster is calculated by incorporating the phase interference of scattered fields from every 
components. Subsequently, the resulting phase matrix is used in the radiative transfer equations 
to evaluate the polarimetric backscattering coefficients from a layer of medium embedded with 
vegetation clusters. Theoretical results are illustrated for various kinds of vegetation clusters. It is 
found that the simulated polarization, frequency, and angular responses carry significant informa- 
tion regarding the structure of vegetation clusters, and also agree with the signatures observed in 
measured multi-frequency polarimetric synthetic aperture radar images. 

In passive remote sensing of earth terrain, the radiative transfer (RT) theory has also been 
widely applied with various types of random medium models to interpret observation. Part of 
the advantages of the RT theory can be attributed to its simplicity in formulating the reflection 
and transmission in the presence of boundaries as well as in treating the scattering effect by the 
incoherent sum of the scattering from each individual scatterer. However, being restricted by this 
incoherent nature, the RT theory cannot properly take into account the coherent effect which 
becomes significant at lower frequency in the case of remote sensing of a structured vegetation 
canopy as observed in our active models. 

In this research, we formulate the vector radiative transfer equation for passive microwave 
remote sensing of a vegetation canopy overlying a soil half-space, and study the calculated bright- 
ness temperatures resulting from microwave thermal emission based on the Gaussian quardrature 
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numerical method. A randomly distributed stem-leave model is employed to construct the phase 
matrix such that the effects of coherence and multi-scale can be properly accounted for. 

For polarimetric remote sensing, geophysical media are modeled as layers containing ran- 
domly embedded scatterers. In a medium such as vegetation canopy, the scatterers can have various 
shapes, sizes, and permittivities that are significantly distinct from the background medium. The 
model studied in this research will consider each type of the scatterers as a species which can takes 
on a shape, size, and complex permittivity different from other species. The effective permittivity 
of the random medium is derived under the strong permittivity fluctuation theory and the po- 
larimetric scattering coefficients are calculated for the layer configuration with the distorted Bom 
approximation in the analytical wave theory which preserves the phase information. 

The multiple species in the random medium are considered as randomly oriented ellipsoids 
and described by multiple three-dimensional ellipsoidal correlation functions. The variances and 
correlation lengths of the correlation functions characterize the fluctuation strengths and the phys- 
ical geometry of the scatterers, respectively. The strong permittivity fluctuation theory is extended 
to account for the multiple species of ellipsoidal shape. In the random medium, a coincidence of 
an observation point with a source point gives rise to the singularity of the dyadic Green’s function 
which is properly taken into account with exclusion volumes of the scatterers. 

Polarimetric scattering properties of a remotely sensed medium are depicted with a co- 
variance matrix whose elements are polarimetric scattering coefficients. The medium has a layer 
configuration: the top layer such as air is considered as homogeneous, the middle layer such as a 
vegetation canopy is random, and the underlying layer such as soil is a homogeneous half space. 
More random medium layers can also be introduced in the configuration to account for weather 
effect such as fog cover. The distorted Born approximation is then used with the effective per- 
mittivity to obtain the polarimetric covariance matrix. The result for the cross-polarized return 
<Tfi V is non-zero even in the first order approximation. Due to the non-spherical shape and the 
random orientation of the scatterers, the correlation coefficient between the HH and VV returns 
has a magnitude differed from unity and a small phase. The scattering coefficients are also used to 
calculate the Mueller matrix for synthesis of polarization signatures. The copolarized signature of 
the random medium has a rather straight distortion track and a recognizable pedestal. 

Polarimetric Passive Remote Sensing of Periodic Surfaces 

To demonstrate the use of microwave polarimetry for passive remote sensing of azimuthally 
asymmetric surface features on the earth terrain, we have performed a series of theoretical analysis 
and experimental measurements for the polarimetric thermal emissions from periodic surfaces. In 
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our recent research, the third Stokes parameter U for thermal emission from periodic surfaces has 
been shown to be zero when the observation direction is parallel or perpendicular to the direction of 
the corrugation, while becomes significant for observation directions between these two directions. 
Thus, the parameter U is related to the relative azimuthal angles of observation and symmetry 
axes of periodic surfaces. This implies that measurement of fully polarimetric emissions of all four 
Stokes parameters can be used to probe the symmetry axes of surface structures. In addition, the 
other measured Stokes parameters are also shown to be related to the height variations of periodic 
surfaces, and permittivities of the underlying medium. 

To verify the theoretical findings, an experiment was designed and carried out by measur- 
ing the polarimetric thermal emissions from a triangularly corrugated soil surface. A triangular 
corrugation pattern was prepared on a surface of soil with known texture and moisture content. 
A 10 GHz radiometer operated at linear polarizations was then used to measure the brightness 
temperatures at different polar and azimuthal angles. The brightness temperatures corresponding 
to the first three Stokes parameters were measured. The fourth Stokes parameter V was not mea- 
sured because the radiometer used could not measure the emission at circular polarizations. The 
measurements were made at horizontal, vertical, and 45° linear polarizations. Significant values 
of the brightness temperature corresponding to the third Stokes parameter U were observed in 
various configurations (as high as 40 K for certain configurations). Theoretical analysis of the mea- 
surement data indicates that the appreciable values of U are caused by the azimuthal asymmetry 
of soil surfaces. 

It is also observed from the experiment that decreases, Tg v increases, Tg p decreases 
to a minimum at <f> = 45° and then increases as takes on the increased values. For Ug y the trend 
is similar to that of Tg p . These general trends are supported by our theoretical predictions of the 
polarimetric brightness temperatures. The significant observation from this experiment is that the 
surface asymmetry can be detected with a measurement of Ug at a single azimuthal angle. 


Composite Volume and Surface Scattering Model 

Among the various theoretical models applied to study the electromagnetic wave scatterings 
from geophysical terrain, such as snow, ice, and vegetation canopy, the radiative transfer theory 
has drawn intensive attention in the microwave remote sensing society during the past years. In 
most of the scattering models, the volume scattering and the surface scattering effects have been 
investigated separately. Recently, there has been a growing interest in the construction of composite 
models which can take into account both types of scattering. 
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In this research, we study the first order iterative solution to the vector radiative transfer 
equations for a two-layer medium with a diffuse top boundary and an irregular bottom boundary 
of Gaussian roughness. The geometrical optics approach with shadowing correction is used in 
formulating the boundary conditions. To demonstrate the utilities of the theory, we apply our 
formula with a phase matrix for randomly oriented spheroidal discrete scatterers to calculate the 
backscattering coefficients from soybean field in different growing stages and compare the results 
with the experimental measurements. Good agreement has been achieved for both the co-polarized 
and the cross-polarized data. It is observed that the presence of the rough surface can significantly 
enhance the backscattering at small incident angles and the levels of the cross-polarized return. The 
polarization signatures calculated based on the Mueller matrix show a straight distortion track and 
an observable pedestal. Numerical comparison to the backscattering coefficients calculated by using 
planar bottom boundary conditions with or without the incoherent addition of the rough surface 
effects are also made in this work. 


Radar Image Classification 

Classification of terrain cover using polarimetric radar is an area of considerable current 
interest and research. A number of methods have been developed to classify ground terrain types 
from fully polarimetric synthetic aperture radar (SAR) images, and these techniques are often 
grouped into supervised and unsupervised approaches. Supervised methods, including both conven- 
tional Maximum Likelihood (ML) and more recent Multi-layer Perceptron classifiers, have yielded 
higher accuracy than unsupervised techniques, but suffer from the need for human interaction to 
predetermine classes and training regions. In contrast, unsupervised methods determine classes 
automatically, but generally show limited ability to accurately divide terrain into natural classes. 
In this research, a new terrain classification technique is introduced, utilizing unsupervised neural 
networks to provide automatic classification, but employing an iterative algorithm which overcomes 
the poor accuracy of other unsupervised techniques. 

Several types of unsupervised neural networks are first applied to the classification of SAR 
images, and the results are compared with those of more conventional unsupervised methods. Neu- 
ral Network approaches include Adaptive Resonance theory (ART), Learning Vector Quantization 
(LVQ), and Kohonen’s self-organizing feature map. Conventional classifiers utilized are the migrat- 
ing means clustering algorithm and the K-means clustering method. With both neural network 
and conventional classifiers, preprocessing is performed to reduce speckle noise and to stabilize the 
training process. Results show that LVQ is the best of the neural network techniques, and that 




this method outperforms all of the conventional unsupervised classifiers. The accuracy of even the 
LVQ technique, however, is seen to remain below that of supervised methods. 

To overcome this poor accuracy, an iterative algorithm is proposed in which the SAR image 
is reclassified using a Maximum Likelihood (ML) classifier. Training of the ML classifier is per- 
formed using a training data set first classified by the above unsupervised method, thus, requiring 
no human intervention, and preserving the unsupervised nature of the overall classification scheme. 
The process is then repeated iteratively, training a second ML classifier using data classified by the 
first. It is shown that this algorithm converges rapidly, and significantly improves classification ac- 
curacy. Performance after convergence is seen to be comparable to that obtained with a supervised 
ML classifier, while maintaining the advantages of an unsupervised technique. 

The new unsupervised and iterative algorithm developed in this research is applied to 
polarimetric SAR images of San Francisco and Beaufort sea ice, acquired by the Jet Propulsion 
Laboratory. The results obtained for this imagery using the new algorithm are compared with the 
results obtained with other techniques, and also with those obtained with single-feature classifica- 
tion. It is found in each case that the new fully polarimetric unsupervised algorithm yields classified 
images which compare closely with those obtained from optimally chosen, supervised algorithms. 
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PUBLICATIONS SUPPORTED BY NASA CONTRACT NAGW-1617 


Refereed Journal Articles and Conference Papers 


Radiative Transfer theory for polarimetric remote sensing of pine forest (H. C. Hsu, H. C. Han, R. 
T. Shin, J. A. Kong, A. Beaudoin, and T. Le Toan,) International Geoscience and Remote Sensing 
Symoposium, Houston, Texas, May 26-29, 1992. 

Application of theoretical models to active and passive remote sensing of saline ice (H. C. Han, J. 
A. Kong, R. T. Shin, S. V. Nghiem, and R. Kwok), International Geoscience and Remote Sensing 
Symoposium, Houston, Texas, May 26-29, 1992. 

Radiative transfer theory for polarimetric remote sensing of pine forest (C. C. Hsu, H. C. Han, R. 
T. Shin, J. A. Kong, A. Beaudoin, and T. Le Toan), submitted to the International Journal of 
Remote Sensing, 1992. 

Scattering of electromagnetic waves from a dense medium consisting of correlated Mie scatter- 
ers with size distributions and applications to dry snow (L. Tsang and J. A. Kong), Journal of 
Electromagnetic Waves Applications, Vol. 6, No. 3, 265—286, 1992. 

Variance of phase fluctuations of waves propagation through a random medium (N. C. Chu, J. A. 
Kong, H. A. Yueh, and S. V. Nghiem), Journal of Electromagnetic Waves Applications, Vol. 6, No. 
2, 169-197, 1992 

Theoretical modeling for passive microwave remote sensing of earth terrain (J. A. Kong), URSI 
Specialists Meeting on Microwave Radiometry and Remote Sensing, Boulder, Colorado, January 
24-26, 1992. 

Polarimetric signatures of a canypy of dielectric cylinders based on first and second order vector 
radiative transfer theroy (L. Tsang, C. H. Chan, J. A. Kong, J. Joesph), Journal of Electromagnetic 
Waves Applications, Vol. 6, No. 1, 19-51, 1992. 

Branching model for polarimetric remote sensing of vegetation (S. H. Yueh, J. A. Kong, J. K. Jao, S. 
Ayasli, and R. T. Shin), Progress in Electromag. Research Symposium, Cambridge, Massachusetts, 
July 1-5, 1991. 

Branching model for vegetation (S. H. Yueh, J. A. Kong, J. K. Jao, and R. T. Shin), Progress in 
Electromag. Research Symposium, Cambridge, Massachusetts, July 1—5, 1991. 

Radiative transfer modelling of vegetation clusters (C. C. Hsu, S. H. Yueh, H. C. Han, R. T. Shin, 
and J. A. Kong), Progress in Electromag. Research Symposium, Cambridge, Massachusetts, July 
1-5, 1991. 
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Application of neural network to polarimetric SAR image classification (Y. Hara, R. G. Atkins, S. 
H. Yueh, R. T. Shin, J. A. Kong, and R. Kwok), Progress in Electromag . Research Symposium , 
Cambridge, Massachusetts, July 1-5, 1991. 

Multiple ellipsoidal species in layer random medium model for polarimetric remote sensing (S. V. 
Nghiem, J. A. Kong, R. T. Shin, and T. Le Toan), Progress in Electromag . Research Symposium , 
Cambridge, Massachusetts, July 1-5, 1991. 

Polarimetric passive remote sensing of an azimuthally asymmetric periodic soil surface (S. V. 
Nghiem, M. E. Veysoglu, J. A. Kong, R. T. Shin, K. O’Neill, and A. Lohanick), Progress in 
Electromag . iZesearcft Symposium , Cambridge, Massachusetts, July 1-5, 1991. 

Polarimetric passive remote sensing of periodic surfaces (M. E. Veysoglu, H. A. Yueh, R. T. Shin, 
and J. A. Kong), Progress in Electromag. Research Symposium, Cambridge, Massachusetts, July 
1-5, 1991. 


Analytical solution of the vector radiative transfer equation with rough surface boundary condi- 
tion (H. C. Hem, J. A. Kong, S. V. Nghiem, and T. Le Toan), Progress in Electromag. Research 
Symposium , Cambridge, Massachusetts, July 1-5, 1991. 

Radiative transfer theory for passive microwave remote sensing (H. C. Han, W. Au, S. H. Yueh, R. 
T. Shin, and J. A. Kong), Progress in Electromag . Research Symposium , Cambridge, Massachusetts, 
July 1-5, 1991. 

Application of random medium model to remote sensing of vegetation (T. Le Toan, S. V. Nghiem, 
J. A. Kong, H. C. Han,) International Geoscience and Remote Sensing Symposium (IGARSS’91), 
Helsinki University of Technology, Espoo, Finland, June 3-6, 1991. 

Application of neural networks to sea ice classification using polarimetric SAR images (R. Kwok, 
Y. Hara, R. G. Atkins, S. H. Yueh, R. T. Shin, and J. A. Kong), International Geoscience and 
Remote Sensing Symposium (IGARSS’91), Helsinki University of Technology, Espoo, Finland, June 
3-6, 1991. 


External calibration of polarimetric radars using point and distributed targets (S. H. Yueh, J. 
A. Kong, and R. T. Shin), Second Annual JPL Airborne Geoscience Workshop, JPL, Pasadena, 
California, May 20-24, 1991. 

Simulation of electromagnetic phenomema using a finite difference- time domain technique (K. Li, M. 
A. Tassoudji, R. T. Shin, and J. A. Kong), 7th Annual Review of Progress in Applied Computational 
Electromagnetics , Naval Postgraduate School, Monterey, 38-55, March 18-22, 1991. 

Electromagnetic scattering of waves by random rough surface: a finite-difference time-domain ap- 
proach (C. H. Chan, S. H. Lou, L. Tsang, and J. A. Kong), Microwave and Optical Technology 
Letters, Vol. 4, No. 9, 355-359, 1991. 

Branching Model for Vegetation (S. H. Yueh, J. A. Kong, J. K. Jao, R. T. Shin, and T. Le Toan), 
IEEE Trans . Geosci . Remote Sensing, Vol. 30, No. 2, 390-402, March 1991. 
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Electromagnetic wave modeling for remote sensing (S. V. Nghiem, J. A. Kong, and T. Le Toan), 
International Conference on Directions in Electromagnetic Wave Modeling^ New York, October 
22-24, 1990. 

Inversion of permittivity and conductivity profiles employing transverse-magnetic polarized mono- 
chromatic data (T. M. Habashy, M. Moldoveanu, and J. A. Kong), SPIE’s 1990 International 
Symposium on Optical and Optoelectronic Applied Science and Engineering, San Diego, California, 
July 8-13, 1990. 

Polarimetric passive remote sensing of periodic surfaces (M. E. Veysoglu, H. A. Yueh, R. T. Shin, 
and J. A. Kong), J . Electromagnetic Waves and Applications , Vol. 5, No. 3, 267-280, 1991. 

Polarimetric passive remote sensing of a periodic soil surface: microwave measurements and analysis 
(S. V. Nghiem, M. E. Veysoglu, J. A. Kong, R. T. Shin, K. O’Neill, and A. W. Lohanick), accepted 
for publication in 7. Electromagnetic Waves and Applications, 1991. 

Application of neural networks to radar image classification (Y. Hara, R. G. Atkins, S. H. Yueh, 
R. T. Shin, and J. A. Kong), submitted for publication in IEEE Trans . Geosci. Remote Sensing, 
1991. 

Calibration of polarimetric radars using in-scene reflectors, (S. H. Yueh, J. A. Kong, and R. T. Shin), 
Progress In Electromagnetic Research , edited by J. A. Kong, Vol. 3, Ch. 9, 451-510, Elsevier, New 
York, 1990. 

Classification and maximum contrast of earth terrain using polarimetric synthetic aperture radar 
images, (J. A. Kong, S. H. Yueh, H. H. Lim, R. T. Shin, and J. J. van Zyl), Progress In Electro- 
magnetic Research , edited by J. A. Kong, Vol. 3, Ch. 6, 327-370, Elsevier, New York, 1990. 

K-distribution and polarimetric terrain radar clutter (S. H. Yueh, J. A. Kong, J. K. Jao, R. T. 
Shin, H. A. Zebker, T. Le Toan, and H. Ottl), Progress In Electromagnetic Research , edited by J. 
A. Kong, Vol. 3, Ch. 4, 237-275, Elsevier, New York, 1990. 

Polarimetric remote sensing of geophysical media with layer random medium model (S. V. Nghiem, 
M. Borgeaud, J. A. Kong, and R. T. Shin), Progress in Electromagnetic Research, edited by J. A. 
Kong, Vol. 3, Ch. 1, 1-73, Elsevier, New York, 1990. 

Application of layered random medium model to polarimetric remote sensing of vegetation (S. 
V. Nghiem, J. A. Kong, and T. LeToan), URSI International Commission F meeting, Hyannis, 
Massachusetts, May 16-18, 1990. 

Calibration of polarimetric radars using in-scene reflectors (S. H. Yueh, J. A. Kong, and R. T. 
Shin), 10th International Geoscience & Remote Sensing Symposium (IGARSS’90), College Park, 
Maryland, May 20-24, 1990. 

Correlation function for a random collection of discrete scatterers (H. H. Lim, S. H. Yueh, R. T. 
Shin, and J. A. Kong), 10th International Geoscience & Remote Sensing Symposium (IGARSS’90), 
College Park, Maryland, May 20-24, 1990. 
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Statistical modelling for polarimetric remote sensing of earth terrain (S. H. Yueh, J. A. Kong, 
R. T. Shin, and H. A. Zebker), 10th International Geoscience & Remote Sensing Symposium 
(IGARSS’90), College Park, Maryland, May 20-24, 1990. 

Study of polarimetric response of sea ice with layered random medium model (S. V. Nghiem, J. A. 
Kong, and R. T. Shin), 10th International Geoscience & Remote Sensing Symposium (IGARSS’90), 
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Abstract- In this paper, the radiative transfer theory is ap- 
plied to interpret polarimetric radar backscatter from pine for- 
est with clustered vegetation structures. To take into account 
the clustered structures with the radiative transfer theory, the 
scattering function of each cluster is calculated by incorporat- 
ing the phase interference of scattered fields from each com- 
ponent. Subsequently, the resulting phase matrix is used in 
the radiative transfer equations to evaluate the polarimetric 
backscattering coefficients from random medium layers embed- 
ded with vegetation clusters. Upon including the multi-scale 
structures, namely, trunks, primary and secondary branches, 
as well as needles, we interpret and simulate the polarimetric 
radar responses from pine forest for different frequencies and 
looking angles. The preliminary results are shown to be in 
good agreement with the measured backscattering coefficients 
at the Landes maritime pine forest during the MAESTRO -1 
experiment. 

' I. Introduction 

In recent years, the application of radar polarimetry for ac- 
tive remote sensing of the earth terrain has inspired extensive 
interests. A variety of theories have been used for electro- 
magnetic modeling of geophysical terrain, such as snow, ice, 
and vegetation canopy [1]. Among these studies, the radia- 
tive transfer (RT) theory has been commonly used to calculate 
radar backscattering coefficients from layered geophysical me- 
dia and to interpret the experimental measurements. 

The radiative transfer theory consists of the radiative 
transfer equations which govern the electromagnetic energy 
propagation through scattering media. Various models have 
been developed based on this theory (2-4). In general, the 
conventional RT theory ignores the relative phase information 
associated with structured scatterers, which may play an 
important role in the overall scattering behavior (5]- In forest, 
vegetation consists of structures of many different scale length 
- the trunk, primary branches, secondary branches, and needle 
leaves, etc. Vegetation elements of each scale are connected 
to elements of other scales in a fashion statistically described 
by the unique architecture pertaining to each tree species. 
For the microwave remote sensing of forest, the vegetation 
*tructures not only give rise to the separation of scattering 
centers for different polarizations, but also provide partially 
coherent scattering by different scatterers with statistically 
prescribed relative positions. 

In this paper, a four-layer RT model is presented for 
the modeling of the pine forest in the Landes area, France^ 
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We make use of the newly developed branching model for 
vegetation [7] to account for the scattering properties of 
structured pine trees. The model thus constructed has a wide 
validity range in frequency spectrum. 


II. Radiative Transfer Theory 


Consider an electromagnetic wave incident upon a multi- 
layered random medium with the incident angle 8 0 as shown 
in Fig. 1 . The scattering regions 1 to 3 are layers of 
thickness h n (n = 1 , 2 , 3 ), where discrete scatterers embedded 
in homogeneous background. The bottom boundary between 
region 4 and region t can be either a fiat surface or a rough 
surface described by a Gaussian random process. 


The vector radiative transfer equation for the specific 
intensity in each scattering region is of the form 


s g WMA = _ ^ (<?| ^ z) . 7 ( 0 , z ) 


J At 


( 1 ) 


where the Stokes vector 7 containing information regarding 
field intensity and phase relation of the two orthogonal polar- 
izations is defined as 


/ = 



1 

V 


( <i^i 2 ) \ 

<|£v| 2 ) 

2Re(E v Ef l ) 

\2lm{E v E' h )J 


( 2 ) 


In (2), the subscripts h and v represent the horizontal and verti- 
cal polarizations, respectively. The angular bracket ( ) denotes 
ensemble average over t he size and orientation distributions 
of scatterers; and 17 = \/a*oAo * ls *^e free-space characteristic 
impedance. 

The extinction matrix ** represents the attenuation due 
to both the scattering and absorption, and can be obtained 
through the optical theorem in terms of forward scattering 
functions. The phase matrix 7(8 i <j>\8 f ,^) characterizes the 
scattering of the Stokes vector from direction into 

(#, <j>) direction. The phase matrix can be formulated in terms 
of scattering functions of the randomly distributed discrete 
scatterers. Along with the boundary conditions, we can solve 
the radiative transfer equations iteratively for the polarimetric 
backscattering coefficients [l]. 
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Figure 2. Cylinder Cluster 


III. Scattering Function for Clustered Structure 

In order to take into account the partially coherent effect 
due to clustered vegetation structures, we formulate the phase 
matrix based on the branching model for vegetation. In pine 
forest, most of the scatterers are of cylindrical shapes. Hence 
in this study, the main subject is cylinder clusters. For the 
cylinder cluster that has one center cylinder and N branching 
cylinders, the total backscattering function / G ^ is 

N 

fa0 = fQap + fnap e (^) 

n=l 

where is the o/?-th element of the scattering matrix for 
the center cylinder, and f nQ g is o/?-th element of the scattering 
matrix for the n-th branching cylinder, a,/?, 7,$ represent 
horizontal or vertical polarization. 

Assuming all the branching cylinders are identical and 
independent of each other, the correlation of / is 

(f<xpfy&) = (foapfoyt) + ^ (fnapfnjS^ 

+ 2N Re [(/oct/?K/n-,|){ e ~*^*)] 

+ N(N - \){f ma p){fn lS ) 

■ <«*'*»)(«-**•) (4) 

The relative phase of the n-th branch with respect to the center 
cylinder is defined as 

— (k t * — fcj) • r n 

where and k s are the incident and scattered wave vectors, 
respectively. r n is the location of the n-th branching cylinder 
relative to the center cylinder. 

In (4), the third and fourth terms are the coherent terms. It 
can be seen that the incoherent approximation is valid when the 
average of the random phase factor (e*^) is so small that the 
coherent terms are negligible as compared with the incoherent 
terms. However, this is not true for vegetation structure with 
scale length comparable with wavelength. 

For the calculation of scattered fields of different cluster 
elements, dielectric cylinders can be used to model trunks, 
branches and coniferous leaves [6,7]. Leaves of deciduous tree 
can be modelled as disks [6,8]. The truncated infinite cylinder 
approximation [9] is employed in this _ study to calculate 
scattered fields from cylinders. 


IV. Comparison with Experimental Data 

The experiment is conducted by the Centre d’Etude Spa- 
tiale des Rayonnemenls (CESR) at the Landes forest, in south- 
west France, using the NASA/JPL Synthetic Aperture Radar 
(SAR) during the MAESTRO-1 campaign in August 1989. A 
strong correlation between P-band backscattering coefficients 
and pine forest parameters has been observed in the measure- 
ment data (10]. 

To model the pine forest, the RT model described in 
previous sections is adopted as depicted in Fig. 1. The top 
layer (region 1) consists of vertically orientated three-scale 
dielectric cylinder clusters charactering the trunks and the 
attached smaller scale branching clusters. These smaller scale 
branching clusters include primary and secondary branches, as 
well as needles. However, at P-band, only the primary and 
secondary branches are important. In the three- scale model, 
all the smaller scale cylinders arc uniformly distributed along 
a larger scale center cylinder with 60° angle relative to the 
orientation of the center cylinder, as show in Fig. 2. Region 
2 contains vertically oriented dielectric cylinders representing 
tree trunks. In region 3, vertical cylinders and randomly 
oriented thin cylinders coexist characterizing a mixture of the 
tree trunks and the short vegetation in the forest understory. 
The underlying grass/ground consists of the attenuating grass 
layer and the ground. 

With this model, we calculate backscattering coefficients 
for six pine stands of different ages. The ages of the forest 
under investigation are 6, 14, 20, 30, 38, and 46 years old. It 
is observed that as forest becomes aged, the average radius 
of tree trunk grows bigger and the average tree height also 
becomes taller. The primary and secondary branches grow 
bigger as well. In general, the above-ground biomass increases 
despite the decrease in the number of trees per unit area 
as forest ages. Since the ground truth indicates that the 
surface is very smooth 'for observation at P-band, we model the 
ground as a planar surface. The grass layer above the ground 
is considered as an homogeneous attenuating layer of 0.2 m 
height and characterized by a dielectric constant (1.05 + i0.5). 
The calculated results for P-band are shown in Fig. 3. The 
discrete points are the data collected by SAR; the curves are 
the theoretical results. It is found from the simulations, the 
main contribution for HH backscattering coefficient is from 
trunk-ground interaction and scattering from branches. As 
for the W and HV, the branches are the dominate scatterers. 
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V. Conclusions 


Radiative transfer theory is applied to the modeling of 
oolarimetric radar backscatter for vegetation canopy. The 
!Lse matrix is derived with the coherent effects due to 
L clustered structures. With this model, we simulate 
multi-frequency and multi-angular backscattering coefficients 
,t different radar polarizations. Theoretical resets are in goo 
agreement with the measurement data collected by NASA/JFb 
SAR- This paper presents the comparison performed at P-Band 
frequency. 
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Abstract- The random medium model is used to interpret the 
polarimetric active and passive measurements of sahne ice. The 
E layer is described as a host ice medium embedded with 
VnHomlv distributed inhomogeneities, and the underlying sea 
u . homogeneous half-space. TTee scauer.m 
11 the ice layer are modeled with an ellipse, dd Ration 
function. The orientation of the scattered is vertically aligned 
S ezimuthally random. The strong permittivity fluctuation 
theory is employed to calculate the effective permittivity and the 
distorted Bom approximation is used to obtain the polarimetric 
Ottering coefficients. We also calculate the thermal emissions 
based on the reciprocity and energy conservation pnnapk^T 
Sets of the random roughness at the air-ice, and ice-water 
interfaces are accounted for by adding the surface scattering to 
the volume scattering return incoherently. The above theoretical 
model which has been successfully applied to analyze the radar 
backscatter data of the first-year sea ice near Point Barrow 
AK, is used to interpret' the measurements performed in the 
CRRELEX program. 

INTRODUCTION 

Understanding electromagnetic wave interaction with saline ice 
prorides physical insight for interpretation of active and passive 
remote senring data of sea ice. Various theoretical models have 
been developed to characterize the scattering and the emission 
properties of sea ice [1]. In this paper, the random medium 
model [2) is used to interpret the measurement data collected 
at the US Armv Cold Regions Research and Engineering Lab- 
oratory (CRREL) under the CRRELEX program. To simulate 
sea ice, saline ice was grown under partially controlled condi- 
tion at CRREL and active and passive microwave signatures 
of the saline ice were measured together with ice characteristic 
parameters. 

Under quiescent condition, the background ice grows in 
columnar foim and saline water is trapped between ice platelets 
in the form of brine inclusions, which are usually ellipsoidal. The 
ice tends to grow vertically downward rendering the ellipsoidal 
inclusions aligned preferably in the vertical direction and t c 
crystallographic c-axos parallel to the Horizontal plane. In t is 
case, the c-axes are, however, random m the horizontal direction. 
The saline ice laver is modeled as a random medium containing 
ellipsoidal scatters vertically aligned but randomly oriented 
in the horizontal direction. The strong fluctuation theory is 
employed to calculate the effective permittivity and the distorte 
Born approximation is used to obtain the polarimetric scattering 
coefficients. The cmissivity of a layer of sea ice overlying sea 
*’fctcr can be obtained through the principles of reciprocity 
and energy conservation [3] In other words, the emissivity 
is obtained as one minus the reflectivity. To account for the 
effects of rough interfaces, we use small perturbation method 
H] to evaluate the surface scattering from air-ice, and ice-water 


interfaces, and add incoherently to the volume scattering return 
obtained with the distorted Bom approximation. This model 
has been successfully applied to study the radar backscatter for 
the thick first-year sea ice data measured near Point Barrow, 
Alaska [5]. In this study, we further employ this model to 
interpret the multi-frequency active and passive measurements 
taken during the CRRELEX experiment. 
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Figure 1 Scattering configuration. 


MODEL DESCRIPTION 

Multi- Layered Random Medium Model 

The fully polarimetric multi-layered random medium model 
for sea icc has been developed with the wave theory. The 
layer configuration is shown in Fig. 1 Region 0 represents 
the homogeneous air above sea ice. The sea ice in region 
contains vertically aligned ellipsoidal scatterers such as saline 
icc with brine pockets or air bubbles which has spatially varying 
permittivity <i(F). Region 2 is the underlying sea water with 
homogeneous permittivity e 2 All three regions are assumed to 
have the free-space permeability n 0 . 

An exponential correlation function locally corresponding to 
a scattcrer is used in the derivation of the effective permittivity 
with the strong fluctuation theory extended to account for 
the ellipsoidal shape and the horizontally random orientation 
of the scatterers. The local correlation function of the brine 
inclusion has three different correlation lengths related to t.ie 
three different axial lengths of the ellipsoid as shown in Fig. 
2. The azimuthal orientation of the scauerer is depicted locaiiy 
with the local z'-axis parallel to the global z-axis and an Eulenan 
fluctuation angle 6/ corresponding to the angle between the c- 
axis and the global horizontal axis (r-axis). When the average 
process is performed over the random orientation angle <?/, 
the effective permittivity is a uniaxial tensor with vertical 
optic axis describing the effective azimuthal symmetry of the 
inhomogeneous sea ice. When a single set of correlation lengths 
is used, for simplicity, to model the brine inclusions with various 
sizes, these correlation lengths are considered as related to the 
effective size of the scatterers. After the effective permittivity 
of the sea ice is calculated, the complete set of polanmetnc 
scattering coefficients axe obtained under the distorted -om 


1265 


31-72SKV92S03.00 © IEE£ 1992 

PRECEDING PAGE BLANK NOT FILLED 


ORIGINAL PAGE IS 
OF POOR QUALITY 




approximation with the anisotropic dyadic Green’s function for 
the layer media. From the scattering coefficients, the Mueller 
matrix is computed to plot the polarization signatures of sea 
icc. 

Under the distorted Bom approximation the polarimetric 
scattering coefficients, obtained as the correlation of the scat- 
tered field, takes the form 


(E 0s {r) £ m 0s {r)) = 

k* j dOj p{oj) J^<&\ C (n j k i m {7\,r^;Of) 

[(c 0 iyp , .n))(rua , i)>] • {(c 0 ut(r,n))(Fi m (n))}' 


(i) 


i 



Figure 2 Geometry of scatterer in region 1. 

where the uniaxia] effective permittivity ! e //l for region 1 is 
used to calculate the mean dyadic Green's functions and the 
mean fields (F\[?))- The correlation function for region 1, 
in (!) is the jT-fm-th element of fourth-rank 

correlation tensor Cc\{ri,T^\<?j) defined as 

C 0jklm (ri,^<Pj) = |W)) P) 

In (1), the probability density function of scatterer orientation 
p(<j>j) is simply 

P(<r>j) “ 7fZ 

for complete random orientation horizontally with no preference 
in azimuthal direction, and ko — Uy/Hoto ^ f ree ' s P ace 
wavenumber with angular frequency u. 


Effective Permittivities 


The strong fluctuation theory is used to derive the effective 
permittivities of the ice layer. The singularities of the dyadic 
Green’s functions in the bilocal approximated Dyson’s equations 
are accounted for; and the low-frequency approximation is 
applied to obtain the results for the anisotropic sea ice media. 
The derivations are done in the frequency domain with the 
Fourier- transform method. 

The ice layer is described with a normalized local correlation 
function of the form 


R iX {r) = exp 


^ ^ Z 1 
7 r *! rJr 7 r 

V V v , 


(4) 


with correlation length tj, and in the local coordinates 
corresponding to the minor, the meridian, and the major axes of 
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the scatterer. In this model, the correlation lengths are related 
to the effective size and shape of the scatterers. The effective 
permittivity of the ice layer is composed of a quasi-static part 
and a scattering part which accounts for the modification in the 
wave speed and attenuation 

? c //i ” ?g\ + ~ W/ 1 ^)] ‘ f«//l (5) 

The effective dyadic scatterer f 4 //i is obtained from the volume 
integration of the dyadic Green’s function and the correlation 
function, the auxiliary permittivity 7 g \ and dyadic coefficient 

5i characterizing the source singularity of the dyadic Green’s 
function are determined by the condition of secular term 
elimination. 


Scattering Coefficients 

For polarimetric backscattering, the scattering coefficients 
are defined by 


O^rt/K — 


lim 

r— dc 
/1-X 


4 ^(E^sEls) 

A Eri EV 


( 6 ) 


where A is the illuminated area, subscript i is for the incident 
wave, s represents the scattered wave, and /i, v, r, and * 
can be h or v which respectively stands for the horizontal or 
vertical components of the electric fields. The components of 
the scattered field in (6) are obtained by measuring the k and 
the v returns while the incident field is transmitted exclusively 
with h or v polarization. For backscattering from reciprocal 
media, the h return obtained by transmitting V is the same as 
the v return obtained by transmitting h. 


Rough Surface Contribution 

To account for the effects of rough interfaces at z = 0 
and z = -d\ t we first approximate the surface statistics by 
a Gaussian correlation function with standard deviation a and 
correlation length £ t then employ the small perturbation method 
to obtain the first-order backscattering field from the interfaces. 
For example, the surface scattering from the air-ice interface is 
given by 


<?hh = ct 2 £ 2 cos 


Ovv = 4fc®CT 2 £ 2 cos 4 6, 
1-2UI-2 


;4 0,|«O 


,TE| 2 -Jfc^sin 2 ^ 


(7a) 


(A:, ~ sin- 9j + fcf.) 


+ klk u ) 7 


e -*|< 5 sin (7b) 


where Oi and Jtj denote the angle of incidence and the wavenum- 
ber for the ice layer, respectively. The surface scattering is inco- 
herently added to the volume scattering return calculated based 
on the method described previously. 


Emissivities 

The Kirchhoff law is invoked to obtain the emissivity of sc* 
ice under isothermal condition from the active bistatic scattering 
coefficients. 


*hV>4>) - i 
*(*,*) = l 


7 - y f 
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where 7bA(0',<^;0,<0) represents the bistatic scattering coeffi- 
cient for an 6-polarized incident wave propagating in the direc- 
tion described by (0,d>) and scattered into the direction 

with 6-polarization. 

DATA INTERPRETATION 

The theoretical results from the random medium model are com- 
pared to the CRRELEX data. For the case under consideration, 
region 0 is air, region 1 is the sea ice layer containing ellipsoidal 
brine inclusions, and region 2 is sea water. The thickness of 
the ice layer is obtained by ground truth measurement. The 
fractional volume of the brine inclusions is calculated from the 
bulk ice salinity and temperature. Results from publications are 
used to estimate the frequency-dependent permittivities of the 
ice background, the brine inclusions, and the underlying sea wa- 
ter. The theoretical calculations performed with and without 
the rough interfaces contributions are compared. 

The local correlation lengths, which are not defined in the 
same manner as the global correlation lengths in the Perovich's 
estimation, corresponding to the effective size of the scatterers 
need to be determined. In our preliminary calculations, instead 
of adjusting the correlation lengths to match the experimental 
data, they are chosen from the characterization of the first-year 
sea ice obtained by a comparison with the 9 GHz measurements 
made at Point Barrow, Alaska [2]. This approach, based on the 
theoretical model, also serves the purpose of referring the remote 
sensing signatures of the saline ice artificially grown at CRREL 
to the sea ice in the natural conditions. 

The theoretically calculated emissivities are in good agree- 
ment with the data. This indicates that the estimated effective 
permittivity tensor, which affects the cmissivity even in the ze- 
roth order, are close to that for natural sea ice At lower fre- 
quencies where the electromagnetic wave can penetrate the ice, 
the emissivities can be sensitive to the ice thickness. In this 
case, the h polarization is more sensitive to the thickness than 
the v polarization which usually has higher attenuation. For 
information pertaining more to the first or higher order effects, 
active signatures should be considered. 

The backscattering coefficients generated theoretically have 
also been compared with that observed in the CRRELEX ex- 
periment for co-polarized and cross-polarized returns only, since 
measurements of other elements in the covariance matrix are not 
available. At 5 GHz. the theoretical results overestimate the CR- 
RELEX data at larger incident angles. This indicates that the 
local correlation lengths of the grown sea ice probably smaller 
than those observed at Point Barrow. At smaller incident an- 
gles. the measurements are higher than the theoretical values 
without rough surface contributions. At 10 GHz. the results, 
emulated using the above correlation lengths, are ai! higher 
‘-••an the measured values provided to us for the backscattering 
* r om the CRRELEX experiment. 
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Scattering of Electromagnetic Waves from a Dense Medium 
Consisting of Correlated Mie Scatterers with Size Distributions 
and Applications to Dry Snow 
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Abstract- The scattering of a plane electromagnetic wave obliquely incident on a layer 
of dense medium consisting of dielectric spherical particles of finite sizes and with size 
distributions is studied. The spherical particles are of sizes comparable to wavelength 
so that Mie scattering is used to describe the single particle scattering characteristics. 
The coherent wave is studied with quasicrystalline approximation using the cross pair 
distribution functions of multiple sizes governed by Percus-Yevick approximation. The 
incoherent scattered wave is calculated with the distorted Born approximation with the 
result expressed in terms of a product of the T-matrices of particles of different sizes and 
permittivities and the Fourier transform of the cross pair distribution functions. The co- 
herent wave effective propagation constants, the attenuation rates and the backscattering 
coefficients are illustrated numerically, with examples chosen to illustrate microwave and 
millimeter wave scattering from snow cover in the frequency range of 5 GHz to 95 GHz, 
and mean grain radius between 0.03 cm to 0.09 cm. Salient features of the numerical 
results for scattering from snow with site distribution are: 1. Correlated dense medium 
scattering is less than independent scattering at low frequency, a fact that is consistent 
with controlled laboratory experiment. 2. Scattering from dense medium of a Rayleigh 
size distribution with an average radius can be much larger than the case of monodisperse 
particles of the sizes identical to that average radius. 3. The scattering attenuation rate 
increases rapidly with frequency at low frequency regime and begins to level off at high 
frequency regime. 4. The coherent wave scattering attenuation rate is large in snow at fre- 
quencies above 15 GHz indicating large optical thickness and the importance of multiple 
scattering. Comparisons with extinction measurements of dry snow at 18 GHz, 35 GHz, 
60 GHz, and 90 GHz are made. 

I. INTRODUCTION 

Propagation and scattering in dense media has applications in remote sensing of 
geophysical terrain and nondestructive evaluation of composite materials [1-13]. 
In a dense medium, particles occupy an appreciable fractional volume. A main fea- 
ture of a medium with densely packed scatterers is that Foldy’s approximation [14] 
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which has been successful in propagation problems through sparse media, is not 
applicable. It has been shown theoretically [1-2,4-6,9,15-16] and through labo- 
ratory experiments [3] that the effects of correlated scattering between scatterers 
must be taken into account. 

Propagation and scattering in dense media has been studied with the quasicrys- 
talline approximation [1-2, 4-6], and the quasicrystalline approximation with co- 
herent potential [6-8] for the first moment of the field and the correlated ladder 
approximation for the second moment of the field [7,9,13]. The dense medium 
theory is developed for media with particles of finite sizes. Because of the finite 
sizes of the particles, their positions axe correlated. The dense medium radiative 
transfer theory has also been developed from these approximations to study mul- 
tiple scattering effects in the incoherent wave. Some of the studies of correlated 
scattering have been limited to the case of low frequency [7,9] and/or particles 
of identical sizes [1-2,4-6,10-11]. However, to approximate physical reality in 
geophysical terrain scattering, size distributions of particles must be taken into 
account. The case of monodisperse particles and poly disperse scattering generally 
can give different results for independent scattering [14,17] as well as correlated 
scattering [8,16]. High frequency results are also important for example at fre- 
quencies above 15 GHz for snow cover. 

In this paper, we study the scattering of a plane electromagnetic wave obliquely 
incident on a layer of dense medium consisting of dielectric spherical particles with 
size distributions. The spherical particles are of sizes comparable to wavelength so 
that Mie scattering is used to describe the single particle scattering characteristics. 
Thus the theory is applicable to particles of arbitrary size. The coherent wave is 
studied with the quasicrystalline approximation using the cross pair distribution 
functions of multiple sizes governed by Percus-Yevick approximation [18-20]. The 
integral equation of the quasicrystalline approximation gives rise to the generalized 
Lorentz-Lorentz law and the generalized Ewald-Oseen extinction theorem [4,21]. 
The generalized Lorentz-Lorentz law consists of a homogeneous system of equa- 
tions. The resulting dispersion relation is an equation governing the effective prop- 
agation constant for the coherent wave. The generalized Ewald-Oseen extinction 
theorem is an inhomogeneous equation that relates the transmitted coherent field 
to the amplitude, polarization, and direction of propagation of the incident field. 
The name “generalized” denotes the fact that the classical Lorentz-Lorentz law 
and the Ewald-Oseen extinction theorem were derived for low frequency effective 
permittivity for the total field with the imaginary part accounting for absorption 
loss only. In the present “generalized” case, the effective propagation constant is 
derived for high frequency when the particle sizes are comparable to wavelength. 
Also the imaginary part of the effective propagation constant includes attenuation 
due to scattering diverted to other directions as well as absorption. 

The calculations of the cross pair distribution functions are based on a spe- 
cial case of Percus-Yevick approximation that has no inter-particle forces except 
for non-interpenetration [8,17-19,22]. The Percus-Yevick cross pair distribution 
functions of multiple sizes have recently been shown to compare very well with 
Monte-Carlo simulations [23]. The incoherent scattered wave is calculated with 
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the distorted Bom approximation with the result expressed in terms of a product 
of the T-matrices of particles of different sizes and permittivities and the Fourier 
transform of the cross pair distribution functions. Numerical examples are chosen 
to illustrate microwave and millimeter wave scattering from snow cover in the 
frequency range of 5 GHz to 95 GHz and mean grain radius ranging from 0.03 cm 
to 0.09 cm [24]. The Rayleigh size distribution and the modified gamma size dis- 
tribution are used to illustrate the scattering. Salient features of the numerical 
results for scattering from snow with Rayleigh size distribution are: 1. Correlated 
dense medium scattering is less than independent scattering at low frequency 
that is consistent with controlled laboratory experiment [3]. 2. Scattering from 
dense medium of a Rayleigh size distribution with an average radius can be much 
larger than the case of monodisperse particles of the sizes identical to that av- 
erage radius. Particles with larger size can contribute to scattering even though 
they are fewer in number. 3. The scattering attenuation rate increases rapidly 
with frequency at low frequency regime and begins to level off at high frequency 
regime. This is to be contrasted with Rayleigh scattering where scattering at- 
tenuation just increases rapidly with frequency. 4. The coherent wave scattering 
attenuation rate can be large at frequencies above 15 GHz based on the Rayleigh 
size distribution indicating large optical thickness and the importance of multiple 
scattering in snow at high frequencies. Comparisons are made with extinction 
measurements of dry snow at 18 GHz, 35 GHz, 60 GHz, and 90 GHz. Good agree- 
ment is obtained and the dense media theory is shown to be able to account for 
the frequency dependence of extinction of dry snow from 18 GHz to 90 GHz. 


II. FORMULATION 

Consider a plane electromagnetic plane wave impinging in direction (n - fa) 
upon a half space of spherical discrete scatterers of multiple sizes and multiple 
permittivities (Fig. 1). Let sj be the species index j — 1,2, ... ,L. Each species is 
described by the permittivity i$j , radius a s j and number density n s j. The num- 
ber density n s j is related to the volume fraction of species f s j, f s j = 4irn s j a V 3 - 
The scatterers are embedded in the lower half space (region 1) with background 
permittivity equal to e which is the same permittivity as that of the region 0 

(Fig. 1). 


The electric field of the incident wave is given by 

E,(r ) = {Ej id + E hi fa)e*« T (1) 

where subscripts v and h denote vertical and horizontal polarizations respec- 
tively, T id denotes the downward propagating incident wave vector and and 
are the associated polarization vectors 

= K sin 0, cos x + sin 0; sin <j>i y — cos 0,z) (2) 

— k%x ^ 1/ ” ^iz ^ 

$id = -x cos 0j cos — y cos 0,- sin <f>{ — z sin 0, (3) 

<f>i = sin <j>i + y cos (f>i (4) 
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In the following, our notation follows that of [6], pp. 446-454 and pp. 506-517. 
Under the quasicrystalline approximation for multiple species, the coherent field 
obeys the equation 

L 


x^ s i)(r 


1 \n ) = 53 n *j f -n) a ( kr i r j)^ 

i ; =l Jv i> 

4. 

i c u mc 


( 5 ) 


where is the conditional average of the exciting field of particles at r/ of 

species s/ conditioning on the particle at rj, T is the matrix representation of 
the T matrix of species sj and <f(Jfcr/Fj) is a transformation matrix transforming 
from vector spherical waves with center at Fj to vector spherical waves with cen- 
ter at F/ as given on p. 452 of [6]. It is an exact wave transformation applicable 
to near field, intermediate field, and far field regions. The function gsjsfrj - n) 
is the cross pair distribution function of two species sj and s/. These can be 
computed readily from the Percus-Yevick approximation. These were originally 
derived for molecular dynamics of fluid mixtures. The result that we use is the 
special case when intermolecular forces are zero, and molecules cannot interpen- 
etrate each other. In (5) the volume of integration Vji for drj is the lower half 
space excluding a “hole” of radius ( a Sj + a Sl ) around vector T\. This is because 

gsjstfj - r,) = 0, for |r; -T { \ < a S] + a s , (6) 



Figure 1. Geometry of the problem: an incident electromagnetic plane 
wave impinging upon a dense medium consisting of spherical 
particles with a size distribution. 
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Scattering of EM Waves from Dense Medium of Mie Scatterers 

In (5) a, nc is the column matrix representing the coefficients of the incident 
wave when expanded into spherical waves, as given on pp. 507 and 514 of [6]. 

To solve (5), we assume the solution 

uK")(r,) = (7) 

where is a column matrix representing coefficients of the exciting field and 

"Kd is the downward going effective wave vector 

K d = K x x + K v y-K z z ( 8 ) 

and the effective wavenumber is 

K = yjK} + K$ + K] (9) 

Substitution of (7) into the integrand of the integral (5), integrations can be 
performed [4,6]. The integral yields two terms that correspond to two waves with 
dependencies 

( 10 ) 

^ e iK t x l +iK y y l -i{k' i -Kl-Kl)' l ! 2 zi (H) 

The first wave given by (10) propagates with K d and is to be balanced with the 
term w s ‘(r t ) on the left hand side of equation since they both have the same 
propagation constant. Balancing these two terms give rise to the generalized 
Lorentz-Lorentz law which can be used to calculate the effective propagation con- 
stant K. The second wave of (11) with incident wave wavenumber k in the square 
root will cancel the second term on the right hand side of (5) which is the incident 
wave term. This cancellation or extinction is the generalized Ewald-Oseen extinc- 
tion theorem [6,25]. For the two terms to cancel out, the wave vector of (11) must 
exactly balance the incident wave vector. This gives rise to the phase matching 
condition 

K x = k iz (12) 

K y = k iy (13) 

so that from (8), (9), (10), we obtain Snell’s law for the effective propagation 
constant, 

K sin0< = k sin 0,- (14) 

K z = (K 2 - k 2 sin 2 8i) l t 2 (15) 

We call both laws generalized because the concept has been generalized to coherent 
wave at high frequency while the classical laws were derived for low frequency case. 

A. Generalized Lorentz-Lorentz Law 

The column matrix for the exciting field can be expressed in terms of 

y.(*>)(W) coe ffi c i en ts where (M) stands for M type of vector 
spherical wave and {N) stands for V type of vector spherical wave and n is the 
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multipole index. The y n ( *> )(M) and y„ ( *> )(N) coefficients satisfy a homogeneous 
set of equations. By setting the determinant equal to zero, one gets a nonlin- 
ear algebraic equation for the effective propagation K. Because the medium is 
statistical isotropic, we have the same homogeneous set of equation for arbitrary 
polarization and incident angle of the incident wave. The homogeneous set of 
equations is 

= -27T £ n„5>n + 1) S p (k t K\R. jt ,) 

*;=1 

{a (1, n| - 1, v\p) A(n, u,p) T^ ,)(M) y„ (M)(,,) 

i(l, n| - 1, v[ p,p - 1) B(n, v,p) Ti' m ^* )( ' ,) } 


+ a| 


(16) 


Y^ W) = — 2tt £ n tl £(2n + l)S p (k, K\R, jtl ) 

*1=1 ”>P 

x {a (1, n| - 1, u\p,p - 1) B(n, u, p) T<" )(M) Y^‘ ){M) 
■ a( l,n| - l, V |pM(n,i/,p)3i' ,)(JV) l r » ( " )W } 


(17) 


where o(l,n| - 1,u\p) “d «(M - l.i'lp.P- 1) are in terms of Wigner 3-j 
symbols and can be found on pp. 449—450 of [6]. In (16), (17) 

v(u + 1) + n(n -f 1) — p(p + 1) 
n (n + 1) 


A(n,v,p) = 


(18) 


. Jin + u + P + l)(y + P ~ n)(n + p - v){n + v ~ P + 1) flW 

*("■ ■'•ri = “ — ^TTj c ; 

Sp{k,K\R $j si) = \khp(kR S} si) jp( K RfjSi) - A hp(kR Sj s t )jp(7 x A* j4( )j 

+ f° dr r 2 [* jfJ (r) - l] h p (kr)j p (Kr) (20) 


— a t] + a Sl 


and 


(*,)(M) _ [p»ijn(p«i)]0~n(C«/) ~ [C«tJn(Cs/)] , Jn(p« < ) 

" “ \p,MP‘i) 1V«(<.,) UG,i»(Cs, )]'*»(*,) 

_ (Pf|jn(P«j)] , C« ) in(Cg<) ~ [Cijjn(C«; )1 pljniPsi) 

” [p«l^n(p*/)] , C 2 1 Jn(C*|) — [C*j,7n(C*J )]' Ps t ^n(Pe t ) 


( 21 ) 

( 22 ) 

(23) 


with p, , = ka f/ , ( '»( = k Sl a t ,, and k t , = , are the Mie scattering T matrix 

coefficients for M and 77 vector spherical waves. 

For each species s/, we terminate at a multipole N s , depending on the radius 
a of that species. Thus (12) and (13) form N c homogeneous equations for 
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the N e unknowns and \ n = 1,2,. . . ,N tl and a/ = \,2,...,L. 

Thus N e = Zf'=iWs r t . . 

The generalized Lorentz-Lorentz law determines the effective propagation K. 
Setting the determinant of (16)-(17) equal to zero gives the nonlinear equation for 
the effective propagation constant K. After K is calculated, then the complex 
transmitted angle 9t and Kz can be determined from (12)— (15). 

The N e coefficients y„ (<,)(M) and Y^ ,){N) are then reduced to only one arbi- 
trary constant that is to be determined by the incident wave. That one equation 
to determine the arbitrary constant is provided by the generalized Ewald-Oseen 
extinction theorem which gives a single inhomogeneous equation. 

B. Generalized Ewald-Oseen Extinction Theorem 

The generalized Ewald-Oseen extinction theorem is obtained by balancing the 
incident wave term of the second term of (5) and the term of (11) that is a result 
of the integral in (5). The result matrix equation can be reduced to one single 

equation for the set of coefficients V„ ( *' )(M) and Y^‘ ){N) . The equation depends 
on whether the incident wave is vertically polarized or horizontally polarized. 


1. Vertical polarized incidence 

In this case E vi ^ 0 and E hi = 0. The coherent transmitted wave is also 
vertically polarized. The generalized Ewald-Oseen extinction theorem gives rise 
to the following inhomogeneous equation 


* ( =1 " 


. (2n + 1) T - ft)) 

|sin(0, — *,)| 


T a; 





(MUj 
2n 1 


(24) 


where the additional superscript ( V ) denotes the vertically polarized case. Note 

that the set of coefficient, and » - 1 JV.„ >, = 

1 L have been determined to within one arbitrary constant by the generalized 
Lorentz-Lorentz law of (16) and (17). Thus (24) provide the last equation that 
determines all those coefficients uniquely. 


2. Horizontal polarized incidence 

For this case E vi = 0 and E hi ^ 0. The coherent transmitted wave is hori- 
zontally polarized. The inhomogeneous equation for the generalized Ewald-Oseen 
extinction theorem is 


. (2rc + 1) 

^ ^ n, **n(n + 1) 
6 { =\ n v 
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+ I " *" V |sin(^i — 0<)| )\ 

_ (K* ~ kj£)h*h. E hi (25) 

where the additional superscript (/?) denotes horizontal polarization. Thus the 
homogeneous system of equations of (16), (17) combined with the inhomogeneous 

equation (25) determines the coefficients and uniquely 

for the use of horizontally polarized incidence. 

The case of an incident wave of arbitrary polarization can be calculated by 
linearly combining the cases of vertical and horizontal polarized incidence. 


C. Bistatic Scattering Coefficients for Incoherent Field 

The incoherent field can be studied with the distorted Born approximation 
or the first order smoothing approximation [4,6]. FYom the incoherent field, the 
bistatic-scattering coefficients for the incoherent field can be calculated. The 
bistatic scattering coefficients are defined by [6] 

t s 4?rr 2 |?Jg| 2 
7*»(*a, *•;*«*) = ( 2^0 cos 0, 


(26) 


where £p is the incoherent scattered field of polarization /?, E l a is the incident 
field of polarization a, j4q is the area of target and r is the distance of observation 
point from target. In the following, we shall give the result for bistatic scattering 
in the plane of incidence <f>& = <f>i 4- ft and 0 < 6$ < tt/ 2. Because of statistical 
isotropy of random distribution of spherical scatterers, bistatic scattering in the 
plane of incidence do not have depolarization for vertically polarized incidence or 
horizontally polarized incidence. Under the distorted Born approximation, 7vv 
and 7 /j/j are given as follows; for 0 < 9$ < tt/2 


I A 1 \ 

7 vv(^«, 4>s = w + <t>i\ 6i, 4>i) = 2^ 2A'"cos0 

* » * 




+ CO s0- n tj n siUtjg t (B t , <j> s — Tt + 4>i\ 4>i) 




WsA) F d 

4ir 

F 


jyr 


(27) 


. . . 4tt Tlsi 

7 hh(9s,<f>s = * + = 72 


2X" cos 0, 




47 r 


: n ») n *l U *)*l( e *' 4>*= z * + <t>u ^i. &) 


A : 2 cos 

X (0 s ,0,)F d 


*1*1 

rW)(t a d.\u/(^)‘f 


(28) 
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lv (V) fn (2n -H) f («,)(M) ( <j )(M)(V) Pn^cosCfl, - <?<)) 

P^(c°s(fl t -g<)) 


T K)(N) y (.,)(^)(V)/ ^ _ P^cos^ - j 

2„ V C ^ * t} |sin(0,-0 ( 


)l 


+ n (n + 1) P n (cos(0, - 0<))) j 

w (W),„ Vr 1 r (2n + *> f T («.)< Ar ) r ^K^)( tf ) P n(cos(g, ~ ft)) 

W « («„«<) = 2J-1) ^Tj\ T " K " | sin(0 4 — 0*)j 

, (^(^(^(MX//)/ _ e) Pn(cos(6s - fit)) 

+ T„ Y n V**"' v t ) | sin( ^_ 0() | 

+ n (n + 1) P„(cos(0 t - 0<))) | 


(29) 


(30) 


and 

U S jsi(0s, 4>s = * + 4>i\ Oil ti) = 2Im(Jf z ) J QQ ^ ~ V c ^ d 




(31) 


where superscript double prime denotes imaginary part. In (31), Kj is given 
by (8), (12)— (15) and 

k s = fc(sin0* cos^*x + sinfl* sin<£ 5 y -fcos0 s z) (32) 

is the wave vector of the scattered wave in region 0. In (31) H $j sp s ji s l = 
1,2,.. , is the Fourier transform of the cross pair distribution function g$jsi 


minus 1. It is also called structure factor. 

too 


= p^s /“ .s'bwcr)- 1]. - **' 


(33) 


Thus the bistatic scattering coefficient is proportional to the Fourier transform 
of the cross pair distribution function, a concept that has been used extensively 
in the investigation of molecular structures of liquids and amorphous solids by 
X-rays [22]. 

For multiple sizes of spherical particles, the cross pair distribution function g S}Sl 
are isotropic, so that SjjSj(F) depends on |r| only. Hence, the cross structure 
factor H S jsi(p) also depends on |p| only. It is useful to note that under the 
Percus-Yevick approximation, closed from expressions are available [18,19,22] for 

HtjSi{p)i &ji s l — 1, 2, • • • , L. 

In (27) and (28) 

F d = (l - e- 7K " d ') (34) 

where K” = Im {K z )- 

We have inserted in (34) the thickness d of the layer containing random discrete 
scatterers. The results of (34) only include volume scattering and ignores any 
possible reflections from the bottom boundary. 
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III. NUMERICAL RESULTS AND DISCUSSION 


In this section, we shall illustrate the numerical results of the effective propagation 
constants and backseat tering coefficients using parameters of snow. The results 
are illustrated with the Rayleigh size distribution which is a special case of the 
modified gamma size distribution. The size distribution is represented by n(a), 
which is the number of particles per unit volume with radii between a and a + da. 


It is given by 


n(a) 


if fa 
16(a) 5 


exp 



(35) 


where (a) 
defined by 


is the mean radius. The fractional volume and the mean radius are 



(36) 


, , J™daan(a) 
{a) ~ f 0 °°dan(a) 


(37) 



Thus, given the fractional volume /, the Rayleigh size distribution contains only 
one parameter, viz., the mean radius. In (38) a m is the mode radius where the 
size distribution is at a maximum (20), and as indicated in (38), is approximately 
equal to 0.798(a). In our numerical calculations of Figs. 2 to 11, we discretize 
the Rayleigh size distribution into 20 different equally spaced sizes (i.e. L = 20 
species). The pair distribution functions for a medium with 20 sizes are then 
calculated by the Percus-Yevick approximation [20]. The maximum size for these 
20 sizes is chosen such that at least 99.7% of the scattering contributions and 
the fractional volume of particles are included. That usually occurs around 2.5 
to 3 times the mean radius. Comparison will be made between the scattering 
of a medium with Rayleigh size distribution of a specified value of (a) and the 
case of monodisperse particles with radii equal to that same value. The results 
of QCA will also be compared with that of independent scattering. For the case 
of independent scattering, the extinction rate is 


2x 


(*e)independent — TJ 


^n 4j 53(2n + l)Re 

.,=1 




(39) 


In Figs. 2 and 3, we plot respectively the effective loss tangent 2K{/ Kt and the 

normalized phase velocity k/Kj- as a function of fractional volume at frequency 

37 GHz for two Rayleigh size distributions with (a) = 0.05 cm and (a) = 0.09 cm. 
We note that the effective loss tangents first increase with fractional volume, rise 
to maxima and then decrease with further increase in fractional volume. These 
are consistent with controlled optical experiments [3] and controlled microwave 
experiments [26]. The peak of the larger size of (a) = 0.09 cm occurs at / = 0.19 
and the peak for (a) = 0.05 cm occurs at / = 0.16. This indicates that scattering 
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by larger particles are “more independent.” “More independent means that for 
scattering by large particles, the “hole correction' for distribution functions very 
often are sufficient to characterize the scattering [6]. Figure 3 shows that the phase 
velocity decreases with fractional volume with smaller velocity for larger particles 
indicating that scattering decreases phase velocity. In Fig. 4, the backscattering 
coefficients are plotted as a function of fractional volume for a layer thickness of 
d = 3 cm. The backscattering coefficients are opp = cosBi with = u,/i. A 
small thickness is chosen because the model of incoherent scattering is limited to 
distorted Born approximation which is valid only for small optical thickness. The 
associated extinction rates of Fig. 2 show that multiple scattering of incoherent 
waves become important for thickness larger than 3 cm as the optical thickness 
becomes large. The results of Fig. 4 indicate that backscattering coefficient can 
be as large as -4 dB for a small layer thickness of d = 3 cm. 



Figure 2. Effective loss tangent 2 Ki/K r of QCA as a function of frac- 
tional volume. Frequency is at 37 GHz. The permittivity of 
the particles is = (3.2 + z0.001)€q. Two cases of Rayleigh 
size distributions are considered with (a) =0.05 cm and (o) — 
0.09 cm. 



Figure 3. Normalized phase velocity k/Kr of QCA as a function of frac- 
tional volume. Frequency is at 37 GHz. The permittivity of 
the particles is €$ = (3.2 4- t0.001)cQ. Two cases of Rayleigh 
size distributions are considered with (a) = 0.05 cm and (a) = 
0.09 cm. 
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Figure 4. Backscattering coefficient Ohh ' n dB of QCA as a function of 
fractional volume, frequency is at 37 GHz. The permittivity 
of the particles is e, = (3.2 + i0.001)e 0 . Two cases of Rayleigh 
size distributions are considered with (a) = 0.05 cm and (a) = 
0.09 cm. The layer thickness is d = 3 cm and incident angle is 

30°. 

In Figs. 5, 6, and 7, we show the coherent wave attenuation rates for a smaller 
mean grain radius of 0.035 cm. The corresponding mode radius is 0.0279 cm. 
Figure 5 compares the extinction rate 2A'; between independent scattering and 
QCA. We note that attenuation rate of independent scattering is larger than QCA. 
The results show that the extinction rate first increases rapidly with frequency as 
in Rayleigh scattering. Then it departs from Rayleigh scattering around 20 GHz 
and starts to gradually level off. This indicates the importance of including Mie 
scattering effects beyond 20 GHz for (a) = 0.035 cm. Figure 6 compares the 
results of QCA between the extinction rates of monodisperse particles with radii - 
0 035 cm and the case of Rayleigh size distribution with (a) = 0.035 cm. We note 
that the case of Rayleigh size distribution has a considerably larger scattering 
rate than monodisperse case. This indicates the importance of including size 
distribution in modelling scattering from realistic geophysical terrain when grain 
sizes are small. It also means that in laboratory controlled experiments when 
particle sizes are made to be close to monodisperse, the measured scattering can 
be substantially less. The numerical results indicate that the extinction relate 
is at 1.056 cm -1 based on Rayleigh size distribution at frequency = 95 GHz 
meaning that there is is strong multiple scattering of incoherent waves at such 
high frequency even for moderate snow layer thicknesses. Figure 7 shows the 
results of HH backscattering coefficient for a layer thickness of d = 10 cm. It 
shows that distorted Bom backscattering can be as large as 0 dB at 60 GHz with 
(a) = 0 035 cm. The curves shows saturation effects due to the limitations of 
distorted Bom approximation. We expect that if multiple scattering of incoherent 
waves are included, the backscattering should continue to increase above OdB. 
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Figure 5. Coherent wave attenuation rate 2 K% in cm"^ as a function 
of frequency. The permittivity of the particles is e s = (3.2 + 
t*0.001)eo. Rayleigh size distribution is considered with (a) = 
0.035 cm and fractional volume / = 0.3. The cases of indepen- 
dent scattering and QCA are compared. 

QCA (single site 1=0.035 ca. vs Rayleigh size <a>=0.035 cm.) 



0 30 40 60 80 100 

frequency in GHs 

Figure 6. Coherent wave attenuation rate 2 A'* in cm” 1 under QCA as a 
function of frequency. The permittivity of the particles is e s = 
(3.2 + i0.001)eQ* The fractional volume is / = 0.3. The cases 
of identical sizes of a = 0.035 cm and Rayleigh size distribution 
with (a) = 0.035 cm are compared. 
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Figure 7. Back scattering coefficient <7/,/, in dB of QCA as a function of fre- 
quency. The permittivity of the particles is e s = (3.2+t0.001)eo* 
The layer thickness is d = 10 cm and incident angle is 30°. The 
cases of identical sizes of g = 0.035 cm and Rayleigh size distri- 
bution with (a) = 0.035 cm are compared. 



In Figs. 8, 9 and 10, the results for the cases of a larger mean grain radius 
of 0.075 cm are shown. The corresponding mode radius is 0.0598 cm. Figure 8 
shows that the extinction rate of independent scattering is still larger than that of 
QCA. However the two curves are much closer than that of Figure 5 showing that 
scattering by larger particles are “more independent.” At frequency = 95 GHz, 
QCA and independent scattering are almost the same. Figure 9 shows that the 
Rayleigh size distribution case is considerably larger than the monodisperse case. 
Strong multiple scattering is to be expected as the extinction rates of QCA for 
Rayleigh size case are large. Figure 10 show the HH backseat tering coefficients of 
a snow layer of thickness d = 1 cm. The small layer thickness is chosen to satisfy 
the validity condition of small optical thickness for distorted Born approxima- 
tion. The saturation effect in Fig. 10 is a demonstration of the limitation of the 
distorted Born approximation. Nevertheless, the backscattering coefficients are 
quite large in view of the small layer thickness of d = 1 cm. In Fig. 11, the VV 
and HH backscattering coefficients are plotted as a function of incident angles at 
frequency = 37 GHz and d = 10 cm. There is a gradual decrease with incident 
angle. There is little difference between W and HH with the HH return slightly 
larger than W. 
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Figure 8. Coherent wave attenuation rate 2K, in cm * as a function 
of frequency. The permittivity of the particles is e s — (3.2 + 
t'0.001)eo- Rayleigh size distribution is considered with (a) = 
0.075 cm and fractional volume / = 0.3. The cases of indepen- 
dent scattering and QCA are compared. 


qca (ilnglc size 1*0.075 cm, aid Kajlelgk »U* <*> * **075 cm.) 



Figure 9. Coherent wave attenuation rate 2K, in cm -1 of QCA as a 
function of frequency. The permittivity of the particles is e s — 
(3.2 + «0.001)eo. The fractional volume is / = 0.3. The cases 
of identical sizes of a = 0.075 cm and Rayleigh size distribution 
with (a) = 0.075 cm are compared. 
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Figure 10. Backscattering coefficient a/,/, in dB of QCA as a function of fre- 
quency. The permittivity of the particles is e s = (3.2+*0.001)eo. 
The layer thickness is d = 1 cm and incident angle is 30°. The 
cases of identical sizes of a = 0.075 cm and Rayleigh size distri- 
bution with (a) = 0.075 cm are compared. 

VV and HH return , d=10 cm., Rayleigh aisc <a>*0.075 cm. 
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Figure 11. Backscattering coefficient a vv and <7/,/, in dB of QCA as a 
function of incident angle. The permittivity of the particles 
is e s = (3.2 + t0.001)f(). The layer thickness is d = 10 cm. 
Frequency is at 37 GHz. The cases of Rayleigh size distribution 
with (a) = 0.075 cm is considered. 


281 


\COCrCC 
\ CILU 

L ou c 
■cl c c < 
Lo< ( c 
c: i c c < 

• o < c c . 
I c c c c < 
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IV. COMPARISON WITH EXTINCTION MEASUREMENTS OF DRY SNOW BE- 
TWEEN 18 GHZ AND 90 GHZ 

In this section, we apply the dense media QCA-PY-Mie theory to compare with 
extinction measurements of dry snow. Dry snow is a mixture of ice grains and 
air. Hallikanen et al. [24] measured the extinction behavior of different types of 
dry snow at four frequencies: 18 GHz, 35 GHz, 60 GHz and 90 GHz. Thus, it is 
important to match the data for a sample over the entire frequency range with 
one set of physical parameters. In this section, we shall match the data of one 
sample (sample number 9 [24]) of large grains. Other samples are being studied 
and comparisons will be reported in the future. 

Recently, grain size distributions of dry snow have been measured extensively 
[27]. Dry snow can be divided into new snow and old snow. New snow has smaller 
mean size and smaller standard deviation while old snow has larger grain size and 
larger standard deviation. The standard deviation of grain size radius can be 
several times larger than the mode radius showing that snow has a broad grain 
size distribution. The maximum grain diameter is between 3 mm to 4 mm. It is 
also observed that the number ratio of large grains to grains with mean size is 
about 1 to 300 [27]. 

In the following, we shall compare with the snow data for sample 9 [24] with 
ground truth of mean grain diameter — 1 mm and snow density = 0.385 gm/cc. 
This corresponds to a fractional volume of 0.423. The reported measurements of 
extinction are 10.2 dB/m, 58.8 dB/m, 247.2 dB/m and 304.0 dB/m respectively 
at frequencies of 18 GHz, 35 GHz, 60 GHz and 90 GHz. Comparisons are made 
with QCA-PY-Mie theory based on two size distributions: the modified gamma 
size distribution and the histogram (Fig. 12) based on information as discussed 
above [27] 

ice grains histogram 
in snow 

1000 1 — 1 



0 12 3 


grain diameter (mm) 

Figure 12, The ice grain size histogram in snow in number per c.c. as a 
function of grain diameter. The histogram is used for QCA- 
PY-Mie theory in Fig. 14. There are 15 sizes, with mean grain 
radius of 0.0479 cm and fractional volume of 0.4134. 
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The modified gamma size distribution is 


n(a) = A'ia p exp (-AT 2 a Q ) 

(39) 


(40) 

* 

ii 

^ CO 

(41) 


Thus, the size distributions have four parameters P, Q, K\, and K 2 . How- 
ever, K 2 and K\ can be expressed in terms of mean radius (a) and fractional 
volume / as indicated in (40)-(41). Thus the alternative four parameters are 

(a),f,P and Q. . 

The second size distribution that we used is in Fig. 12 which shows the his- 
togram number density (number per c.c.) for 15 sizes with mean grain radius 
of 0.0479 cm and fractional volume of 0.4134. Note that there are large grains 
between 3 mm to 4 mm with the number ratio to grains of mean size of about 1 
to 300 [27]. 

In Fig. 13, we make a comparison with the data based on the modified gamma 
size distribution with P = 1,Q = l A, (a) = 0.05 cm. and / =0.423. The permit- 
tivity of ice grains chosen are e* = 3.2«o + *c* with e”/eo = 0-007, 0.009, 0.011 
and 0.014 respectively at 18 GHz, 35 GHz, 60 GHz and 90 GHz which corresponds 
to a salinity of 0.12 parts per thousand [28]. The results of independent scattering 
are also shown for comparison. The results in figure 13 show that the independent 
scattering results are larger than the data while the QCA-PY-Mie theory shows 
a reasonable match based on the modified gamma size distribution. 

In Fig. 14, we make a comparison with data based on the histogram size distri- 
bution of Fig. 12 and using the same permittivity of ice grains as in Fig. 13. The 
agreement is good. It seems that with size distribution, the QCA-PY-Mie the- 
ory can explain the frequency dependence of the extinction of dry snow between 
18 GHz and 90 GHz. The results of independent scattering are also shown for 
comparison. Future comparisons will be made when the size distrbutions of ice 
grains in snow become available [27]. Note that we have used the Percus-Yevick 
pair functions for non-interacting particles without mutual force. It seems that 
ice grains may have sticking force that may affect the pair distribution functions. 
Besides QCA, the quasicrystalline approximation with coherent potential [6,8] 
should be used. Monte-Carlo simulations of extinction rate can also be made 
based on the solution of Maxwell’s equations [29]. 


Scattering of EM Waves from Dense Medium of Mie Scatierers 


283 




c • cU c 

C <_ V. w 

C 9 . 

C< 

0 9 C 

CH 

t, 9 : 


Figure 13. Comparison between QCA-PY-Mie theory and independent 
scattering with snow extinction data at 18 GHz, 35 GHz, 
60 GHz, and 90 GHz of Hallikanen et al. [24] with ground truth 
of mean grain diameter of 1 mm and snow density =0.385 gm/ cc. 
The theory is based on the modified gamma size distribution 
with P = 1, Q = 1.4, < a >= 0.05 cm, and / = 0.423. 
The permittivity of ice grains chosen are e = 3.2 cq + ie* with 
o = 0.007,0.009,0.011 and 0.014 respectively at 18 GHz, 
35 GHz, 60 GHz, and 90 GHz which corresponds to a salinity of 
0.12 parts per thousand [28]. 

ground truth of snow: 1mm mean 
grain size; 0.385 gm/cc 



Figure 14. Comparison between QCA-PY-Mie theory and independent 
scattering with snow extinction data at 18 GHz, 35 GHz, 
60 GHz, and 90 GHz of Hallikanen et al. [24] with ground truth 
of mean grain diameter of 1 mm and snow density=0.385 gm/cc. 
The theory is based on the histogram of Fig. 12 with < a >= 
0.0479 cm, and / = 0.4134. The permittivity of ice grains cho- 
sen are — 3.2 cq -I- ie$ with c^/e o = 0.007,0.009,0.011, and 
0.014, respectively, at 18 GHz, 35 GHz, 60 GHz, and 90 GHz 
which corresponds to a salinity of 0.12 parts per thousand [28]. 
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V. CONCLUSIONS 

In this paper, we have illustrated the importance of dense medium theory of 
correlated Mie scattering for studying snow cover at microwave and millimeter 
wave frequencies above 15 GHz. We have also demonstrated the importance of 
including size distributions as there are large difference between the polydisperse 
and the monodisperse cases. The importance of polydispersions has long been 
noted [17], particularly at low frequency when scattering is proportional to particle 
size to the sixth power. Based on size distributions, the dense medium theory is 
able to account for the frequency dependence of extinction of dry snow between 
18 GHz and 90 GHz. The calculations also show that extinction rates are large 
for dense medium so that multiple scattering effects become important even for 
small layer thicknesses. The calculation of the incoherent wave in this paper has 
been limited to that of distorted Born approximation. The subject of multiple 
scattering for the incoherent wave including the effects of correlated Mie scattering 
should be an important future subject for snow cover scattering beyond 15 GHz. 
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THEORETICAL MODELING FOR PASSIVE MICROWAVE REMOTE SENSING OF 

EARTH TERRAIN 
by 
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Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139, USA 
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Abstract 

In this paper we address the issue of theoretical modeling for passive 
microwave remote sensing of terrain media such as snow, ice, vegetation, 
and periodic surfaces. Historically, the volume scattering effects 
stimulated the development of the continuous random medium model and 
the random discrete scattering model for the description of the media. 
Theoretical treatments were developed along two different paths. 
Invoking the principle of reciprocity, the wave theory based on Maxwell's 
equations has been used to calculate the emissivity. The other approach 
was to start with the radiative transfer equations and solved for the 
brightness temperatures directly. Attempts have been made to derive the 
radiative transfer theory from the wave theory. At the same time, both 
theoretical approaches have been used to calculate the radiometric 
emissions and to interpret experimentally measured data. 

The successful interpretation of the Cosmos 243 data was perhaps the 
first most important step towards a serious development of the 
continuous random medium model to account for the volume scattering 
effects of snow ice fields. Subse’quent interpretation of measurement 
results from snow field with both passive radiometers and active radar 
systems established a unique position for its description of earth terrain 
media. Recent efforts in classifying sea ice with correlation function 
characterization are demonstration of acceptance -of this model. Future 
inverse scattering developments will perhaps rely heavily on this model. 



In this paper, we shall illustrate the development of the theoretical 
models and present data matching results with measurements made in 
snow ice fields and vegetation canopies. The emissivity calculations for 
periodical rough surfaces will also be presented and compared to 

measured data. 

Recent development in polarimetric active remote sensing with synthetic 
aperture radar has created significant theoretical results and ^pracHcal 
applications. In passive remote sensing, the third and the fourth Stokes 
parameters for earth remote sensing have not received much attention in 
the past partly due to the expected small values any measurement can 
yield We have made initial calculations and experimental measurement 
to show that at least the third Stokes parameter can give appreciable 
number in both theoretical prediction and actual experimentation in a 
plowed field. Such results may have practical implications in measuring 
wind directions in ocean waves, for instance. We believe polarimetric 
passive remote sensing is a viable field which should be explored in light 
of its potential applications that may derived from the full Stokes vector 
instead of its first two parameters, i.e.,ters, i.e., the horizontally and 
vertically polarized brightness temperature components. 

Acknowledgement: The work presented in this paper contains 

contributions from L. Tsang. R. Shin, S. Nghiem. H. Yueh, J. Lee. M. Borgeaud. 
K. O’Neill. M. Veysoglu, H. Han. C. Hsu, and others. 
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Im/ uv is given by (20) below. In the second case, consider the incident wave to 
be of unit amplitude in the horizontal polarization, i.e., E v = 0 and E \ = 1. 
We then calculate the interior field of the particle i?£(r / ) and then the imaginary 
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